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Abstract 

Three-dimensional random tessellations that are stable under iteration (STIT 
tessellations) are considered. They arise as a result of subsequent cell division, 
which implies that their cells are not face-to-face. The edges of the cell-dividing 
polygons are the so-called I-segments of the tessellation. The main result is an ex- 
plicit formula for the distribution of the number of vertices in the relative interior 
of the typical I-segment. On the way of its proof other distributional identities for 
the typical as well as for the length-weighted typical I-segment are obtained. They 
provide new insight into the spatio-temporal construction process. 

Keywords. Cell division process; iteration/nesting; marked point process; ran- 
dom tessellation; stability under iteration; stochastic geometry 
MSG. Primary 60D05; Secondary 60G55, 60E05. 

1 Introduction 

In recent years, random tessellation theory has been an active field of research. Whereas 
in the past mainly mean values and their relations were considered, current research 
focuses on second-order parameters, limit theorems and distributional results, see [U 
[5l |6] to mention only a few. Besides the classical Poisson hyperplane and Poisson- 
Voronoi tessellations, random tessellations constructed by subsequent cell division have 
attracted particular interest in recent times in stochastic geometry and spatial statistics 
(see [miKH] and especially [3] and the references cited therein). Among these models, 
the so-called STIT tessellations (which are stable under iteration, see below) introduced 
in [9l|T2] are of particular interest, because of the number of analytically available results 
[UlinillS], IIS]-[IS] and L2QJ-[23J, see Figure □ for illustrations. The model shows the 
potential to become a new reference model for crack or fissure structures. 

After a detailed analysis of planar STIT tessellations in [TUl [ISl ED] , in the present pa- 
per the 3-dimensional case is considered. We study so-called I-segments in homogeneous 
three-dimensional STIT tessellations, which form the one-dimensional building blocks of 
a STIT tessellation. They appear in the course of the sequential cell splitting procedure 
when cells are divided by new planes. In fact, all sides (1-faces) of a dividing two- 
dimensional polygon - which will be referred to as an I-polygon - are called I-segments. 
In principle the distribution of the number of vertices in the relative interior of the typ- 
ical I-segment is in the focus of this paper. Regarding the spatio-temporal construction, 
it becomes evident that this number depends on the birth time, the direction and on 
the length of the I-segment. But it also depends on the birth time of that I-polygon, in 
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Figure 1: Realisations of an isotropic STIT tessellation (left) and of a STIT tessellation 
whose directional distribution TZ is concentrated with equal weight on three orthogonal 
directions (right). 



whose interior the I-segment under consideration arises at a later time. This polygon 
is called the carrying I-polygon of the I-segment. At its birth time, this I-polygon is a 
facet of two adjacent mosaic cells and during the cell division process these cells can 
undergo further subdivision and in this way new vertices and edges within the carrying 
I-polygon can appear, see Figure [21 

For this reason we consider the edges in the carrying I-polygon (dark grey) and 
observe that an I-segment can have intersections with already existing edges (light col- 
oured lines in the dark grey polygon) at the moment of its birth, see Figure [2j This is 
in sharp contrast to the planar case [W\ [2U] and causes that their analysis is consider- 
ably more involved. Further vertices can arise after the birth of the typical I-segment 
as an effect of the ongoing cell division procedure. For deriving an explicit formula 
for the probability that a fixed number of vertices is located in the relative interior of 
the typical I-segment, we need to study at first the marked process of I-segments with 
the following marks: length, direction and birth time of the I-segment and birth time 
of the carrying I-polygon. We will describe the joint and all marginal distributions of 
the marks associated to the typical and for the length-weighted typical I-segment of a 
homogeneous spatial STIT tessellation. Moreover, by noting that such a spatial tessel- 
lation has exactly two different types of vertices (see Figure ED we can refine our result 
and obtain the joint distribution of the numbers of vertices of both types in the relative 
interior of the typical I-segment. 

Note that it would also be interesting to calculate the distribution of the number of 
vertices within the typical 2-dimensional building block, i.e. the typical I-polygon. How- 
ever, this quantity is currently not accessible, because to this end the area distribution 
of a typical Poisson polygon is needed. To determine this distribution is a long-standing 
open problem in stochastic geometry. 

As already mentioned, related problems in the planar case have been studied in 
[iniEH- The results obtained there have led to a deeper understanding of planar STIT 
tessellations and were used to obtain new results about their structure. In a companion 
paper [22] we also use the results of the present work to explore in more detail the 
combinatorial structure of spatial STIT tessellations, a study initiated in |2H I23j. 

The paper is organized as follows: In Section [2] we introduce some notation, recall 
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Figure 2: An I-segment (bold) together with its carrying I-polygon (dark grey) at the 
time of its birth (left) and after further subdivision (right). 



the basic construction of STIT tessellations and rephrase some of their most important 
properties, which are needed for our later arguments. The conceptual framework as well 
as our main results are contained in Section [3l The proofs are the content of Section HI 

2 STIT tessellations in 

In this section we explain our basic notation, the construction and the main properties 
of STIT tessellations that are needed for our arguments below. Although our results are 
dealing with the three-dimensional case, we focus in this section on general space dimen- 
sions, because in our proofs we apply some of the properties also to lower-dimensional 
tessellations. 

2.1 Notation 

A tessellation of with d > 1 can be described in two ways: as a locally finite 
collection of non-overlapping and space-filling compact convex polytopes (called cells in 
the sequel), as well as a closed subset of M"^, which is formed by the union of all cell 
boundaries. We switch between both perspectives arbitrarily. 

In this paper we will deal with random tessellations of W^, which can be regarded as 
random variables taking values in the measurable space of tessellations of W^, see jlStllOj 
for a detailed definition and in particular for measur ability issues. Our attention will 
be restricted to homogeneous (spatially stationary) random tessellations, i.e. random 
tessellation whose distribution is invariant under spatial translations. 

We denote by T-L the set of all hyperplanes in M'^. A hyperplane h ^ % will be 
parametrized by its normal direction u € Sj^^ (the upper unit half-sphere in W^) and 
its signed distance p € M to the origin, where the distance is defined as positive iff the 
projection of the origin o onto h is located in the upper half-space. Such a hyperplane 
is denoted hy h = h{p, u) € Ti. For i? C M'^ define 

[B] = {{p, u) eRx S^-^ : h{p, u) n S ^ 0} 

as the set of parameter values (p, u) of hyperplanes h{p, u) hitting the set B. Consider 
a measure A on M x Sj^^ (equipped with the Borel product cr-field), which is the image 



3 



under the described parametrization of a (non-zero) locally finite, translation invariant 
measure on H. Invariance under translations implies that A factorizes, i.e. there is a 
constant A > and a probability measure TZ on Sf'^ with 

(1) A = A//®7^, 

where /i denotes the Lebesgue measure on R. If TZ is the uniform distribution on 5^^^ 
and A = 1, then A corresponds to the isometry-invariant hyperplane measure Aiso, cf. 
|15j . From now on we assume A = 1 in the factorization Moreover, to avoid 

degenerated cases we assume that TZ is not concentrated on a great half-subsphere of 
i.e. we require span(support(7^)) = M"'. 

2.2 Construction 

A formal and detailed description of STIT tessellations in bounded windows is given 
in [12j . Here we explain the construction in an intuitive way only. Let A be a hyper- 
plane measure as above and W C M.'^ a bounded convex polytope. We assign to a 
random lifetime and on expiry of this lifetime, we choose a random hyperplane, which 
splits W into two polyhedral sub-cells and W~. The construction continues now 
independently and recursively in both of the sub-cells and W~, which is to say 
that and W~ are provided with independent random lifetimes and that they are 
divided by random hyperplanes when they die. Note that the hyperplanes are always 
chopped-off by the boundary of their respective mother cells. This repeated cell division 
is continued until a fixed deterministic time threshold t > is reached. The random 
tessellation constructed until time t within W is denoted by Y{t, W). 

In order to ensure the temporal Markov property of the described cell splitting 
process, we assume that the lifetimes of the cells are conditionally (given the cells at a 
certain time) independent and exponentially distributed. Moreover, we assume that the 
parameter of this exponential lifetime distribution of a cell c is given by A([c]) (and thus 
the parameters of the exponential distributions of different cells are not independent). In 
the special case A = Aiso we have that Aiso([c]) is proportional to the integral-geometric 
mean width of c. Note that this choice ensures that smaller cells live stochastically 
longer. Furthermore, we will assume that the hyperplane splitting a cell c is chosen 
according to the law A([c])~-^A(- n [c]). 

2.3 Important properties 

We summarize here those properties of STIT tessellations that are needed in our argu- 
ments below. For further background on STIT tessellations we refer the reader to [9l [T2] 
and to [15l [19] for a general introduction to stochastic geometry. 

If we insert independent copies Yc{s) of the tessellation Y{s) into the cells c of Y{t) 
(generating Yc{s) fl c), the resulting iterated or nested tessellation is denoted by 

Y{t)mY{s) = Y{t)U \J {Y,{s)nc). 

c a cell of Y{t) 

Considered as a random closed set, Y{t, W) satisfies the following 

[Spatial consistency] The random tessellation Y(t, W) is spatially consistent in that 

for any convex V C W with W as in Subsection [ZD we have Y{t, W)nV = Y{t, V) 

(where = stands for equality in distribution). Thus, Kolmogorov's extension the- 
orem ensures that there exists a homogeneous random tessellation Y{t) in the 

whole satisfying Y{t) nW = Y{t, W). For the particular choice A = Aiso, Y{t) 
is also isotropic, i.e. its distribution is rotation invariant. 
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In [9] a global construction of Y(t) is provided, where the spatio-temporal random 
process (y(t))t>o is defined. Its important properties are summarized here: 

[Scaling] The distributions of the rescaled tessellations are identical, i.e. tY{t) = 
sY{s) for all s,t>0. 

[Iteration stability] We have 

(2) Yit) = Y{s)mY{t- s), 0<s<t, 

where Y[s) and Y{t — s) are independent. Consequently, because of [Scaling], 

Y{t) = nY{nt) = n{Y{t) ffl . . . ffl Y{t)) 

' V ' 

n times 

for t > and n € N. The latter relation is usually referred to as stability un- 
der iterations. For this reason, the random tessellations Y{t) are called STIT 
tessellations. 

[Poisson typical cell] The distribution of the interior of the typical cell of Y(t) 
coincides with that of the interior of the typical cell of a homogeneous Poisson 
hyperplane tessellation with intensity measure tA (cf . [151 [T9] ) . 

[Interpretation of t and TZ] The surface density Sy of Y{t), this is the mean total 
{d — l)-volume of cell boundaries of Y{t) per unit d- volume, is equal to the con- 
struction time t > 0, i.e. Sy = t. The probability measure TZ, see ([1]), is the 
distribution of the normal direction in the typical boundary point of Y{t), called 
the directional distribution of the tessellation (this is the surface-area-weighted 
directional distribution of the cell boundaries). 

[Linear sections] The intersection of Y{t) with a line L induces a homogeneous Pois- 
son point process on L with intensity A([e(L)])t, where e{L) is a segment of unit 
length on L. 

[Temporal Markov property] The random process (y(t))t>o with values in the 
space of tessellations satisfies the Markov property in time. 

3 Framework and results 

In this section we present our main results in Theorem [H Theorem [5l Theorem [7] and 
their corollaries. Before we can state them, we have to introduce further notation and 
the framework of our work. 

3.1 The marked process of I-segments in spatial STIT tessellations 

The notion of I-segments is due to R. Miles and characterizes one of the types of line 
segments associated with a random planar tessellation that is not face-to-face, cf. 
and the references cited therein. An I-segment can be defined as a maximal union of 
collinear and connected line segments that appear in the edge-skeleton of a tessellation. 
Obviously, this definition does not depend on the dimension of the ambient space and 
will henceforth also be used for the spatial case, d = 3. Generalizing Miles' concept, 
an I-polygon is a maximal union of coplanar and convex planar polygons (subsets of a 
plane in R^). Here, 'maximal' is understood in the sense that there is no extension to 
a larger collinear (coplanar) and convex set in the tessellation (considered as the closed 
set of cell boundaries). 
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In terms of the spatio-temporal construction of STIT tessellations Y{t), the I- 
segments and I-polygons in spatial STIT tessellations can be characterized as follows. 
Any I-polygon of Y{t) is a cell-dividing planar polygon introduced until time t, and 
any I-segment is a side (1-face) of an I-polygon. When a cell c is divided by a plane 
h € [c], exactly one I-polygon cDh is born, whereas at least three I-segments are born 
simultaneously. Any of these I-segments is the intersection of the I-polygon cDh with a 
facet (2-face) of the cell c. This facet itself is embedded in a (possibly larger) previously 
born I-polygon, which is called the carrying I-polygon of the I-segment. 

For our purposes it will be appropriate to describe the joint distribution of the 
following marks associated to an I-segment in Y{t): 

(3) {£, ^, P, /?carr) G (0, Oo) X 5^ X (0, t) X (0, t). 

Here i, (p and (3 are the length, direction and the birth time of the I-segment, respectively, 
and /Scarr is the birth time of the carrying I-polygon. Note that the direction of a segment 
is the unique unit vector in parallel to the segment. 

Consider for fixed time t > and measure A the homogeneous STIT tessellation 
Y{t). Then the process of the marked I-segments of Y{t) with marks as in ([3]) is a 
homogeneous marked segment process. Thus, Palm calculus for marked point processes 
can be applied and allows to define the distribution of the typical I-segment and its 
mark distribution, see [151 [19] a-^d also Section H] below. In intuitive terms, the typical I- 
segment and its mark distribution can be regarded as a randomly (equally likely) chosen 
I-segment of Y{t, W) together with its random marks when is a 'large' observation 
window. 

In the following we use the indicator function notation 1{. . .}, which is 1 if the 
statement in brackets is fulfilled and otherwise. Let us also define the two constants 

(4) C2 := I l[u,v]n{du)n{dv), 

si si 

u,v,w] TZ{du)Tl{dv)TZ{dw) , 




si si Si 



where [u, v] is the area of the parallelogram spanned hy u,v € Sf_ and where [u, v, w] 
stands for the volume of the parallelepiped spanned by u,v,w G 5^ (interpreted as 
unit vectors in having one endpoint at the origin). Note that in the isotropic case 
these constants are given by C2 = f and Cs = f • We further introduce the following 
probability distributions on S^: 

(5) TZ{u) ■■= ^ [ [ i{u^nv^nsleu}[u,v]n{du)n{dv), 



si si 



(2 

(6) 7^typ(c/) - § / e u}A{[u]) n{du), 

'si 

where U C is a Borel set and where A([u]) stands for the A- measure of the set of 
all planes hitting the unit line segment connecting u € 5^ C M'^ with the origin o, and 
U"*" denotes the orthogonal complement of u € . For a Poisson plane tessellation with 
intensity measure A the distributions TZ and T^typ are the length-weighted distribution 
of the direction of edges and the directional distribution of the typical edge, respectively, 
see In the particular isotropic case we have A([ii]) = ^ for any u G 5^ and thus, TZ 
and T^typ are the uniform distribution on S'^. 



6 



3.2 Statement of results 

We are now prepared to present the main results of this paper. We start with a descrip- 
tion of the distribution of marks, see Q, associated with the typical I-segment: 

Theorem 1. Let Y{t) be a homogeneous random STIT tessellation in M'^ with measure 
A as in 

(i) The distribution of the direction of the typical I-segment of Y{t) equals T^typ 

defined in 

(ii) The joint birth time density of (/3, /3carr) of the typical I-segment of Y{t) equals 

Pp,Pc.rAs,r) = ^1{0 < r < s < t}. 

(iii) The conditional length density of the typical I-segment of the STIT tessellation 
Y{t), given (v?, /3, /3carr) = {u,s,r) G 5| x (0, t)^ with < r < s < t, is 

P£|^=«,/3=.,/3_=r(^) = A(M)^ e-'^(M)-l{x > 0}. 

It is interesting to note that the conditional length density in part (iii) of the previous 
theorem does not depend on r and the joint distribution of the birth time vector (/3, /3carr) 
in part (ii) is independent of A (resp. TZ). 

Now we derive some consequences of Theorem [TJ At first, we calculate some marginal 
and conditional birth time densities: 

Corollary 2. For the homogeneous random STIT tessellation Y{t) we have 

PiSc.rri'r) = -^—l{0<r<t}, 
(7) p[s{s) = — l{0<s<t}. 

Moreover, 

2s 

P/3|/3,,„=r(s) = ^1{0 < r < s < t}, 

P/3_|/3=,(r) = ^l{0<r<s}, 

i.e. the conditional birth time distribution of the carrying I-polygon, given that the I- 
segment is born at time < s <t, is the uniform distribution on (0,s). 

Using now Theorem [T] (ii) and (iii) we obtain 

Corollary 3. The direction (f and the birth time pair (/3,/3carr) of the typical I-segment 
of Y{t) are independent. Moreover, the joint conditional density of length and birth 
times, given its direction ti G 5^ of the segment, is 

P£,A/3_|^=n(^,^,0 = 3^(M)"' e-A([n])sx > q^^^q < r < S < t}. 

Integration yields now the marginal density for the length of the typical I-segment 
of the STIT tessellation Y{t). This density is already known from [18], see also Lemma 
[T3] below: 
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Corollary 4. The length density of the typical I-segment ofY(t) equals 



+3A{[u])H'^x^ + A(['u])3t3a;3)g-A([«])te^ TZtypidu), x > 0. 
In the isotropic case this reduces to 
, . 3 



Pi[x) 



t^x^ 



48 - (48 + 2Atx + Qx^t^ + t^x^)e 



X > 0. 



Now we turn to the distribution of the number of vertices in the relative interior of 
the typical I-segment. In fact, with Theorem [1] we obtain the following: 

Theorem 5. The probability p„ that the typical I-segment of Y[t) has exactly n G N 
vertices in its relative interior is given by 



(8) 



1 1 

P„=3//(l 





(3-(i-.)(3-t,)r ^^^^ 

^ (3-(l-a)(2-6))"+i 



It is interesting to note that the distribution in Theorem [5] does not depend on the 
measure A (or equivalently the directional distribution TZ) and the time parameter t. 
But, this is evident, because the number of vertices in the relative interior of the typical 
I-segment does not change when the tessellation is scaled in space. However, the latter 
is, because of [Scaling] , equivalent to a rescaling of time. Some particular values for p^, 
are summarized in the following table. 





Pn (exact value) 


p„ (numerical value) 


n 


= 


i|2ln3-261n2-f 


0.43289 


n 


= 1 


ifl^ In 3- 107 In 2 -35 


0.21384 


n 


= 2 


^In3- 3501n2-^ 


0.11841 


n 


= 3 


31617 In 3 1025 In 2 ^^99 


0.07075 



From formula ^ in Theorem [5] we conclude 

Corollary 6. The mean number of vertices in the relative interior of the typical I- 
segment is equal to 2. Moreover, the variance equals 59/3 and all higher moments of 
that random variable are infinite. 

Note that the mean value is in accordance with the result in |21] . Furthermore, the 
non-existence of higher moments is not surprising, because also moments of order > 3 
of the length of the typical I-segment of Y{t) are infinite, too. 

As a last result we would like to point out that Theorem [5] admits a refinement, 
which is needed in [22] . We note at first that a spatial STIT tessellation has exactly two 
diff'erent types of vertices called T-vertices and X- vertices, see |21)-|23). Illustrations of 
these two types of vertices are shown in Figure [3l Given a carrying I-polygon (dark 
grey), a T- vertex is generated if two further I-polygons intersect in the same half-space 
determined by the carrying I-polygon. An X-vertex is generated by an intersection of 
two further polygons in the two different half-spaces specified by the carrying 1-polygon. 

For the typical I-segment we denote by Pm,n, m, n € N, the probability that it has 
exactly m vertices of type T and n vertices of type X in its relative interior. 



8 



Figure 3: A T-vertex (left) and an X-vertex (right) in a spatial STIT tessellation. 



Theorem 7. For a homogeneous spatial STIT tessellation and m, n G N we have 

-3.2-f"^+"'i / 1(1 A- a-a-«)a-fe)r .^^^ 

Pm,n-6 I y ^ jj j{l a) a (3_(-^_^)(2-6))-+"+i'^'^ 



From this formula the following can be concluded: 

Corollary 8. The mean number of T-vertices and of X-vertices in the relative interior 
of the typical I- segment equals 1, respectively. Furthermore, the variance of the number 
of T-type vertices is 8 and that of the number of X-type vertices is 11/3 (all higher 
moments are infinite). In addition, the covariance of the number of T- and X-vertices 
equals 8. 

In |2H [23] it is shown that the proportion of the intensities of T- and X-vertices in a 
spatial STIT tessellation is 2 : 1. Whereas any T-vertex is located in the relative interior 
of exactly one I-segment an X-vertex is in the relative interior of two I-segments. This 
confirms the first statement of Corollary [51 



4 Proofs 

Before proving our results for the typical I-segment, we first consider the mark distri- 
bution in a typical edge point, i.e. the length-weighted mark distribution. The corres- 
ponding results are derived in Subsection 14.21 The proofs of Theorem [H Theorem [5] and 
Theorem [7] and their corollaries are the content of Subsection 14.31 Some preparatorial 
material is collected in Subsection 14.11 



4.1 Some preparations 

We denote by £ the measurable space of line segments in W^, and equip it with the Borel 
cj-field induced by the Hausdorff distance, cf. ^15j. The distribution on £ of the typical 
I-segment of a STIT tessellation Y(t) for fixed t > and fixed measure A is denoted by 
D^*^*^ and, similarly, the distribution of the typical edge of a Poisson plane tessellation 
with intensity measure sA will be denoted by B>^^'^\ 
In Theorem 3 in [18] the following is shown: 

Lemma 9. For any non-negative measurable function f : C M we have 

t 

I f{L) B^('\dL) = ^f{L) D^W(dL)ds. 
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In order to apply this proposition we will need the length distribution of the edges in 
Poisson plane tessellations. We consider the survival function of the conditional length 
distribution of the typical edge of P{s), given its direction n € 5^, i.e. 



H, 



rP{s) 

£\(p=u 





use a tilde to indicate that we refer to a length-weighted distribution. 

Lemma 10. The conditional length distribution of the typical edge in the Poisson plane 
tessellation P(s) has the survival function 



and the survival function of the conditional length distribution of the length-weighted 
typical edge in the Poisson plane tessellation P{s) equals 



Proof. The intersection of the homogeneous Poisson plane process P{s) with a fixed 
plane parallel to u € 5^ is a Poisson line process. The intersection of this line process 
with a fixed line parallel to u is a Poisson point process with intensity A([ii])s, see e.g. 
[15\ Theorem 4.4.6] or, more specifically, (4.31) ibidem. Thus a twofold application of 
Slivnyak's theorem for Poisson processes yields that the typical edge point of a Poisson 
plane tessellation is a.s. located on a length-weighted segment that is generated by a 
linear homogeneous Poisson point process, and, under the condition that the direction 
is u, its intensity is A([m])s. Hence the length of the typical edge, under the condition 
ip = u is exponentially distributed with parameter A([ii])s, and the respective length 
weighted distribution is the Gamma distribution with the parameter (2, A([ii])s).n □ 

As a corollary we obtain the mean length of the typical edge in a Poisson plane 
tessellation. 

Lemma 11. The mean length of the typical edge in the Poisson plane tessellation P{s) 



Proof. With Theorem 1 in [6j , which says that T^typ is the directional distribution of the 
typical edge in a Poisson plane tessellation, and ([6]) we obtain 




for X > 



<il(:E) = (l + A(M)5x)e-^(M)^^- for x > 0. 



IS 



fP(s) _ 1^ 



oo 



oo 







Hp''{x)dx 




£\ip=u 



{x)TZtyp{du) dx 



oo 




si 



which proves our claim. □ □ 



We derive now from Lemma a similar representation for the distribution B^^*) of 
the length-weighted typical I-segment of Y{t) in terms of D^^''^, the distribution of the 
length-weighted typical edge of a Poisson plane tessellation. 
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Lemma 12. For any non-negative measurable function f : C ^ R we have 

t 

I f{L) m'\dL) = ^f{L) m'\dL)ds. 



Proof. Denoting the mean length of the typical I-segment of Y{t) by ^* and that of the 
typical edge of a Poisson plane tessellation P{s) by we have 

j f{L)W^\dL) = i |/(L)length(L) B^W(dL) 

t 

= 1/ |^/(L)length(L)D^(^)(dL)ds 



t 

= i l^-fF(^^^lf{L)m^\dL)ys 



with Lemma [9l Using now ^* = from |21] and Lemma [TTl we conclude the assertion 
of the lemma. □ □ 

Now we consider the survival function of the conditional length distribution of the 
typical I-segment of Y{t), given its direction u € i.e. 

Hli^=ui^) = J l{length(L) > x} ©^(^^^^JdL) 

Y(t) 

where 1} ^^j^^^^{dL) denotes the respective conditional distribution of the typical I- 

segment. Analogously, the length-weighted case H^i\^^^ is defined. The two distribu- 
tional identities in Lemma[9]and Lemma [12] together with Lemma [TO] implv the following 
result. 

Lemma 13. The length distribution of the typical I-segment ofY{t), given its direction 
u € S\, has the survival function 

t 



and the length distribution of the length-weighted typical I-segment of Y{t), given its 
direction u E 5^, has the survival function 

t 

H\\^=u{^) = / f (1 + ^^{[u])sx)e-^^^^^>^'ds. 



4.2 Length-weighted mark distributions 

For construction times < r < s we consider the states Y{r) and Y[s) and introduce 
three different random sets of unions of I-segments of Y{s). These are: 

(a) the union of all I-segments of Y{s), 

(b) the union of all those I-segments of Y{s) with birth time in (r, s\ that appear in 
the interior of the cells of Y{r), 
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(c) the union of I-segments of Y{s) with birth time in (r, s] that appear on the facets 
of the cells of Y{r). 

For these random sets the corresponding length measures, the length intensities and 
the directional distributions in the typical edge point (the length-weighted case) are, 
respectively, denoted by 

(a) ^(s, .), Lv{s) = E^(s, [0, and %, 

(b) - r, .), L*y{s -r)= Efi*{s - r, [0, l]^) and Q*/"', 

(c) ij{r, s, .), Lv{r, s) = E^l{r, s, [0, l]^) and Q'/ . 

For a Borel set B CM.^, e.g. ^*(s — r, B) is the total length of edges of Y{s) with birth 
time in (r, s], which are located in the interior of the cells of Y(r) and in B. It is evident 
that, for < r < s, we have for the length measures 

(9) //(s, •) = /i(r, •) + //(r, s, •) + /i*(s - r, •)• 

From the STIT property the following relations can be deduced: 

Lemma 14. We have Ly{s — r) = Lv{s — r) and Ly(s) = ^2*^ with C,2 given by 
Furthermore, Q*/''' = Q*"*^ = Q^^ = TZ with U defined via 



Proof. Resorting to <^ we find Ly{s — r) = Ly{s — r) and Q^'*~' = Q^"''- Equation 
Ly(s) = C,2S^ is (8) in [13] and the fact that does not depend on s follows from [131 
Eq. (14)]. The remaining equality = 7^ is a consequence of Eq. (11) ibidem. □ □ 



In a next step we consider for the STIT tessellation Y{t) with t > the joint length- 
weighted distribution ^ ^^^^^ of direction G 5^ and birth times /?, /3carr £ (0, t]- 

Lemma 15. For any Borel set U C Sl_ and < r < s < t we have 



QUp^^JU X {r,s] X (0,r]) = ^ilf^n{U). 

The corresponding joint density of (/3,/3carr) with respect to the Lebesgue measure on 
{0,tf is 

_ 2 

Pl3A.Js,r) = ^1{0 <r<s<t}, 

i.e. (/? and (/3, /3carr) o-f^ independent and (/3, /3carr) 'is uniformly distributed on the triangle 
{(s,r) G M2 : < r < s < t}. 

Proof. We use ([9]), to conclude 

(10) Ly(s) =Ly(r)+Lt.(s-r) + Ly(r,s) 

for the length intensities defined above. In view of Lemma Owe obtain from (jlOp . 

(11) Ly(r,s) =2C2(rs-r2). 

From the definition of the length-weighted directional distribution it follows that the 
mean total length per unit volume of I-segments of Y{t) with direction in a Borel set 
U (ZS\ is Lv{t)n{U). Combining this with ([TUD and ^ yields 

Ly(r, s)'R{U) = 2C2(rs - r^)Tl{U). 
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Thus, for the mark distribution in a typical edge point we find 

QU,^_ju X (,.,,1 X (o,.i) = "-'2;g"'' = '-^m 

and hence, ip and (/3,/3carr) are independent. Since /3carr < /3 with probability 1 we have 
for the joint distribution function Ff^^0^^^^{s,r) of /3 and /3carr with < r < s < t the 
equation 

(12) q'^^f,_{{r,s] X (0,r]) =F^,^_(s,r)-Q*,,^_((0,r] x (0,r]), 

where /3^^j,j,((?', s] x (0, r]) = ^ /3carr('^+ ^ ^ '^1) stands for joint distribution 
of (/3,/3carr) in a typical edge point of Y{t). Partial differentiation with respect to r 
and then with respect to s yields Pi3,i3c^^^{s,r) (note that the Q-term on the right hand 
side of (fT2]l does not depend on s and thus vanishes after differentiation w.r.t. s) and 
completes the proof of the lemma. □ □ 

We turn now to the joint length- weighted distribution ^ ^ ^^^^^ of length, direction 
and birth times of the I-segment, i.e. of the I-segment through a typical edge point of 
Y{t). The key is the conditional length distribution of the I-segment, given its direction 
u € 5^. For any point z in the edge-skeleton of y(t) denote by £{z),ip{z), /3{z), PcsLvriz) 
the a.s. uniquely determined mark of the I-segment through z. According to the defin- 
ition of the mark distribution (see p. 84 in |15]) we have for any Borel set U C 5^ using 
the Campbell theorem and ([9]) 

Ql^,/3Aa..((^'°o) X [/X (r,s] X (0,r]) 
^ E / l{i{z) > x,(p{z) G U}fi{r,s,dz) 



Lvit) 
1 

MO 



E 



l{f!.{z) > x,ip{z) G U}fi{s,dz) - J l{£{z) > x,(f{z) G U}fi{r,dz) 

l{i{z) > X, Lp{z) G U}fi*{s - r, dz) 
^^^^ ^Lv{s)qi^{{x,<^) X U)-Lvir)qi^iix,^) x U) 

(13) -Ly(s-r)Q*;^-'^((x,oo) X [/) 

where Q^'^"'' is the {£, (/9)-marginal distribution of Q*e'^~i3 i3^^^^- Write now Q£|^=„ for the 
length-weighted conditional distribution of length under the condition ip = u. Then 

Qljix, oo)xU) = y 1{/ > X, n G U}qi^{d{l, u)) 

(14) = J ll{l>xW,^^^^{dl)l{ueU}nidu), 

where we have used Lemma UM In view of (|13|) and p4|) we now have to determine 
certain conditional length distributions. 

Lemma 16. For 0<r<s,0<x and a Borel set U C we have 
Lv{s)Ql^{{x, oc) X U)-Lv{r)Ql^{{x,oo) x U) 
= C2 [ A( J)2^2 ([3 + Mlu])rx(3 + A([n])rx)]e-^(M)- 



u 



A{[u])^x^ 

[3 + A([u])sx(3 + A(M)sx)]e-^(M)^^) TZidu). 
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Proof. We start by recalling from LemmafT^that Lv{s) = C,2S^ and Ly{r) = C,2i''^- Using 
now ()14p and the integral representation for the conditional survival function provided 
in Lemma [13] we calculate 

Lv{sm,^{{x,^) X U)-Lv{rni^{{x,^) x U) 
u u 

s 

= 2C2J y"r;(l + A([M])ra)e-^(M)™d?;^((iu) 

U r 

and integration proves the claim. □ □ 

It remains to determine the last item Q^'^~^((x, cxo) x U) in ()13p . which in view of 
Lemma [13] may be written in the form 

(15) Qlf-'-{{x,oo) xU) = j jl{l> x}Q*';i;(d/)l{n e U}R{du). 

Lemma 17. For 0<r<s, 0<x and a Borel set U C we have 

2 



Lv{s - r)Q^;^ {{x, oo)xU) = C2j ^^^jp^a 

u 

X (3(e-^(["l)™ - e-^(["])"^) + A([n])2sx2(re-^(M)™ - se-^(M)^^)) n{du). 

Proof. To derive the result we make use of a method developed for the planar case in 
[8] and consider the conditional survival function of the length of the remaining 

I-segment, this is the part of the I-segment that lies above a random point of the edge 
skeleton, which is selected according to ^*(s— r, •). (In this proof, the letter G will always 
refer to a remaining I-segment while H refers to the whole I-segment.) Regarding 
the distribution of remaining length is the distribution of the minimum of the remaining 
length in Y{s — r) and the distance (in direction u) to the 'frame' tessellation Y(r), 
where the I-segment is cut. Because of [Linear sections] the latter has the conditional 
survival function ^-Hi'^iy^ _ Making use of the independence of Y{r) and Y{s — r) we 
obtain 

(16) G;|';rjx) = G,^|-;_(x) • e-^(N)-, 

where G\~^-^ is the corresponding conditional survival function for the remaining I- 
segment of the tessellation Y[s — r). Palm theory yields (compare with Eq. (14) in 
[8]) that the survival function of the length of the remaining I-segment as considered 
above and the conditional survival function H^^^_^ of the length of the (whole) typical 
I-segment are related by 

00 

(17) G^r" (x) = -V f H'r' {a)da, 



where i^^u is the conditional mean length of the typical I-segment in Y(s — r). Using 

V=« ~ 2A([uJ)(s-r) • 



Lemma [13] we calculate ^2;=^ = 2MU IVs-r'i • Combining ([16]) with ([T7]) and using again 



Lemma [13] we find 



00 s—r 



(18) G^jrjx) = ^^e-^(M)- j j v^^^^l-^y-dvda. 



X 
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We use now once more the Palm theory developed in [8j (see in particular Eq. (16) 



there) to conclude that the survival function Hl'^J'^ of conditional length distribution 



of the corresponding typical I-segment equals 

and thus we find 

= A(M)»(»-"p(2» + r).3 + A(|„l)rx)e-<1"1)" 

-(2 + (2s - r)A{[u])x + K{Mf<s - r)x2)e-^(M)^^ 



from p8p by integration. Using this formula we calculate the corresponding mean I- 



segment length C^Lu as 



n*,s — r / zj*,s-'r I 



A([n])(2s + r) 



and again by length-weighting we get the corresponding survival function 



oo 



which is given by 

T*,s—r / 



A(M)2sa;2(re-^(M)'^" - se'^^M)^")) . 



Taking into account the factor Lyi^s — r) = (^2(3 — r)"^, we can complete the proof by 
using ([I5|). □ □ 



Combining now ()13p with Lemma [T6| Lemma [T7I and the fact Lv{t) = C2t'^ from 
Lemma [13] we arrive at 

Ql^,/3Aa..((^'°o) xUxir,s]x (0,r]) 

2r 



u 



((2 + A([n])rx)e-^(["1)''^ - (2 + A([n])sx)e-^(["1)^^') 7^(d^z). 



A relation similar to ()12p for ^ ^ ^^^^^ and the respective distribution function yields 
after differentiation w.r.t. r, s and x the following assertion: 

Lemma 18. The conditional joint density of length i and the birth time vector (/3,/3carr) 
of the length-weighted typical I-segment in Y{t), given its direction = u £ S'^, equals 

i^.,/3Aa.|^=n(^> r) = Mi[!^:,e-^(M)- l{x > 0}1{0 <r<s<t}. 

We turn now to the proofs of our main results. In fact, Lemma [18] is the key for the 
proof of Theorem [T] which forms the basis to derive Theorems [5] and [7] 
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4.3 Proof of the main results 

Proof of TheoremUl (i) Apply Lemma [9] with f{L) = \{'^{L) € C/} for a measurable 
U C 5^ and where '^{L) the direction of L. Since for homogeneous Poisson plane 
tessellations with intensity measure sA, s > 0, the directional distribution of the typical 
edge is invariant w.r.t. the scaling factor s, this implies that the distribution of the 
direction of the typical I-segment of Y{t) is the same as the directional distribution of 
the typical edge in a Poisson plane tessellation with intensity measure tA. Thus we can 
apply Theorem 1 in [6], which immediately yields assertion (i). 
(a) The key is the relation 

^*=«~ 

where is the conditional mean length of the typical I-segment in Y{t), given ip = u. 
With Lemma [13] we calculate = 2tK{[u]) ' ^'^'^ Lemma [18] implies 

(19) = ^A(H)se-^(["l)^^ > 0}. 
Integration with respect to x yields now 

oo 

(20) P/3,/3_|^=«(^r) = ^1 A([u]).e-^(M)-da; = ^, < r < . < t. 



Since it does not depend on u, we have pp^p^^^^ = Pi3,fs^^„\(p=u ^iid this completes the 
argument. 

(in) Similarly to the proof of part (ii) above, we use the following relation for the 
conditional densities: 

ft 

Pe\^=u,l3=s,l3c^rr=r[X) - Pe\^=u,l3=s,Pc.rr=r[X), 

where ^^^^^^^^ is the conditional mean length of the typical I-segment in Y{t), 

given if, (3 and /3carr- Combining (fT9]l and ([20|) we calculate 

Pei,^u,,.s,,....Ux) = ^^f^-"'^^"^";;';^ = A(M).e-A(N)^n{. > 0}, 

which implies ^^^^ ^^^^^^j. = A(Jn[)s' This completes the argument. □ □ 

Proof of Corollaries The formulas in Corollary [2] follow by straightforward integ- 
ration from Theorem [1] (ii) . 

To see Corollary [3] we first note that in the proof of Theorem [1] (ii) we have already 
shown that 

P.,/3Aa.|^=n(x,s,r) = -^A(M).e-^(M)-, 

which proves the formula. Integration with respect to x yields now the conditional joint 
density of the birth time vector, namely 

( \ _ 
Pl3A^rAv=u{s.,r) - 

This shows the independence of if and the birth times (/3,/3carr) and finally completes 
the proof of Corollary [3] 

The statement of Corollary H] follows by integration from the formula p9|). □ □ 
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Proof of Theorem O The first step is to verify the following expression for 

t s oo 

^H[u])s^ sMlu])x 
t3 




(21) 52 



n\ 

We start by noting that an I-segment arises as the intersection of two I-polygons, which 
are born at two different time instants, the maximum of which is the birth time of the 
I-segment and has a density pp{s) given by d?]). However, in contrast to the planar 
case, I-segments in a spatial STIT tessellation can already have vertices in their relative 
interior at their time of birth. These vertices are generated by extant I-segments 'on the 
backside' of the carrying polygon, see Figure [2] (left). Given its length ^ = x G (0, 00), its 
direction (/9 = u G 5^, its time of birth /3 = s S (0, t) and the birth time /3carr = r € (0, s) 
of the carrying I-polygon, the number A'^birth of vertices in the relative interior of the 
typical I-segment at the time of birth s has, by ^ and [Linear sections], a Poisson 
distribution with parameter A([n])x(s — r). Moreover, in the time interval (s,t], on 
both sides of the I-polygon with birth time s and 'behind' the carrying I-polygon a 
Poisson distributed number of vertices appear in the relative interior of the segment, 
whose parameter is, by ([2]) and again by [Linear sections], A([u])x(f — s), see Figure 
[2] (right). Thus, the sum of these independent numbers is again Poisson distributed, 
but with parameter 3A([ii])x(t — s). Adding the independent number A^'birth leads to a 
Poisson distribution with parameter A([n])3;(3t — 2s — r). 

Recall now that we have calculated in Corollary [3] the conditional joint density of 
length and birth time vector and, thus, mixing the Poisson distributed number of points 
with respect to this density and the directional distribution T^typ, we arrive at ([2T]l . 

We show now that (|2ip is equivalent to ([8]). First, using 

00 

(A(M)x)"e-^(M)^"e-^(M)"(3t-2s-r)^^ ^! 





we see that (1211) can be transformed into 



A([M])(3i-s-r)"+i' 



t s 

3s2 (3t-2s-r)" 



52 

* ' 3s2 (3t - 2s - r)" , , 

aras. 



+1 



{3t- s- r)" 


Applying now in the inner integral the substitution 6 = 1 — r/s, we obtain 

/ /"Ss^ (3t-s(3-6)r 

''" = 7 7 l^ (3t-s(2-b))^+i ^^^^ 


and applying the similar substitution a = 1 — s/t, finally leads to ([8]). □ □ 
Proof of Corollary {Bl The numbers shown in Corollary [6] can immediately be calculated 
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from the explicit formula in Theorem [5j To see it, write 



1 1 



Vnp, = 3/ hl-afYn }^ ~ ~ ~ dbda 
^ ^ 11^ ' ^ (3 - (1 - a)(2 - 



n=0 Q Q n=0 

1 1 



3 y" J{l-a){a{3-b) + b)dbda = 2 







and similarly 



2 } }r ^3V^ 2 (3-(l-a)(3-6))" „ , 



n=0 g Q n=0 

1 1 



(3- (l-a)(2-6))"+i 
y j a^{2b'^ - lib + 15) - 0(46^ - 106 - 3) + 6(1 + 2b)da db 




71 

which yields the variance value 71/3 — 2^ = 59/3. □ □ 

Proof of Theorem Let the birth time /3 = s < t of the typical I-segment and the birth 
time /3carr = r < s of the carrying I-polygon be given. At time s, the typical I-segment 
can only have X-type vertices in its relative interior (this is the number A^birth in the 
proof of Theorem [5] above) and further vertices of this type can be created from 'behind' 
the segment until time t. At time s the I-segment is an edge (1-facet) of two cells. 
Vertices of type T in the interior of the I-segment can appear in the time interval (s,t) 
by further division of these two adjacent cells. After these observations the proof readily 
follows the lines of the proof of Theorem [5] and for this reason the details are omitted. 
□ □ 

Proof of Corollary O Corollary [8] follows from Theorem [7] by direct calculation. For 
example, the mean number of T-type vertices in the relative interior of the typical 
I-segment can be calculated as follows: 



00 00 00 I, x3 



m=0 n=0 m=0 



m-O g g 

1 1 

6 J a{l — a)da J db = I. 



dbda 



For the second moment we find 



1 1 



m=0 71=0 m=0 

1 1 

= 6 y a(l + 3a)(ia J db = 9 


so that the variance equals 9 — 1^ = 8. Similar calculations yield mean 1 and variance 
11/3 for the number of X- vertices in the relative interior of the typical I-segment. The 
value for the covariance is obtained from Corollary [6] by calculating 59/3 — 11/3 — 8 = 8, 
which completes the argument. □ □ 
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